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The large majority of neutron stars (NS) in low mass X-ray binaries (LMXBs) have never shown 
detectable pulsations despite several decades of intense monitoring. The reason for this remains an 
unsolved problem that hampers our ability to measure the spin frequency of most accreting NSs. 
The accreting millisecond X-ray pulsar (AMXP) HETE J1900. 1-2455 is an intermittent pulsar that 
exhibited pulsations at about 377 Hz for the first 2 months and then turned in a non-pulsating source. 
Understanding why this happened might help to understand why most LMXBs do not pulsate. We 
present a 7 year long coherent timing analysis of data taken with the Rossi X-ray Timing Explorer 
(RXTE). We discover new sporadic pulsations that are detected on a baseline of about 2.5 years. We 
find that the pulse phases anti-correlate with the X-ray flux as previously discovered in other AMXPs. 
We place stringent upper limits of 0.05% rms on the pulsed fraction when pulsations are not detected 
and identify an enigmatic pulse phase drift of ~ 180° in coincidence with the first disappearance of 
pulsations. Thanks to the new pulsations we measure a long term spin frequency derivative whose 
strength decays exponentially with time. We interpret this phenomenon as evidence of magnetic field 
burial. 

Subject headings: pulsars: individual (HETE J1900. 1-2455) — X-rays: binaries 



1. INTRODUCTION 

Some neutron stars in LMXBs have magnetic fields 
which are sufficiently strong to truncate the accretion 
disk and channel plasma along the field lines. The NS ro- 
tation modulates the X-ray emission emerging from the 
hot spots plus the accretion shocks that form close to 
the NS surface. Detecting their spin has several impor- 
tant implications for understanding how millisecond ra- 
dio pulsars are recycled, how the magnetospherc and the 
accretion disk interact and whether sub-millisecond NSs 
can be formed. 

Accreting millisecond X-ray pulsars in LMXBs spin 
with periods of less than 10 ms and are powered by chan- 
neled accretion. Only a small fraction of neutron stars in 
LMXBs are AMXPs (see Patruno 2010b for the AMXP 
list and Papitto et al. 2011 for the most recent system 
discovered), the largest majority not showing accretion 
powered pulsations with a fractional rms amplitude (see 
Eq. 4 in Hartman et al. 2008 for a definition) smaller 
than ~ 1% (Vaughan et al. 1994; Dib et al. 2005; Patruno 
2010a). The reason for this is still a puzzle and several 
models attempt to provide an explanation by invoking 
gravitational lensing(Ozel 2009), pulse smearing in a hot 
electron cloud (Titarchuk et al. 2002), rotation and mag- 
netic pole alignment (Lamb et al. 2009a, b), MHD insta- 
bilities (Romanova et al. 2008) and burial of the magnetic 
field (Cumming et al. 2001). 

This paradigm remained almost unchanged until 2007- 
2008, when three AMXPs showed intermittent pulsa- 
tions: HETE J1900. 1-2455 (Kaaret et al. 2006; Galloway 
et al. 2007), Aql X-l (Casella et al. 2008) and SAX 
J1749.8-2021 (Gavriil et al. 2007; Altamirano et al. 2008; 
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Patruno et al. 2009a). In these systems X-ray pulses ap- 
pear and disappear on timescales that range from a few 
hundred seconds up to several days. In this respect the 
discovery of intermittent AMXPs has been an important 
breakthrough since it was suggested that most of the non- 
pulsating LMXBs might indeed sporadically pulsate. 

Particularly interesting is the behavior of HETE 
J1900. 1-2455 which showed persistent pulsations for the 
first 22 days until the occurrence of a flare in the 
lightcurve on July 8, 2005 (MJD 53,559). Right after this 
event the pulsations disappeared and reappeared at dif- 
ferent intervals with the source now becoming intermit- 
tent until August 20, 2005 (MJD 53,602). After this date 
pulsations disappeared and they were only tentatively de- 
tected with fractional amplitude of 0.29% rms, by sum- 
ming the power spectra of 137 ks of data (Galloway et al. 
2008). This is so far the smallest fractional amplitude 
ever reported for an AMXP. This slow transition from a 
normal AMXP to a non-pulsating LMXB might there- 
fore be the key to understand why most LMXB do not 
pulsate. 

The outburst of HETE J1900.1-2455 has lasted for - 7 
years since the first discovery on June 15, 2005 and is still 
ongoing at the moment of writing this Letter. We present 
a coherent timing analysis of the 7 years of data collected. 
We identify an enigmatic 180° drift in the pulse phases 
discovery during the July flare, and we identify new pul- 
sations in a few data segments that extend to a few years 
after the last robust detection on August 20, 2005. We 
use these new pulsations to measure the behavior of the 
spin frequency derivative over a baseline of 2.5 years. We 
then discuss how these findings might help to understand 
the pulse formation mechanism and the behavior of non- 
pulsating LMXBs. 

2. X-RAY OBSERVATIONS AND DATA ANALYSIS 
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We used all high time resolution data taken during the 
lifetime of the RXTE with the Proportional Counter Ar- 
ray. We used data modes with time resolution of 2~ 20 
s (GoodXenon) and 2 -13 s (Events-122/^ s). We selected 
an energy band which spans approximately the range 2- 
16 kcV (absolute channels 5-37) to maximize the signal- 
to- noise ratio (S/N) of the pulsations. The data are 
barycentered using the JPL DE405 ephemeris at the best 
determined optical position of HETE J1900.1-2455 (Fox 
2005) and are cleaned according to standard procedures 
with X-ray bursts removed from our analysis. 

Pulsations are constructed by folding the data in seg- 
ments of length 300, 500, « 3000 seconds (i.e., orbit-long 
RXTE observations), a few to several hours (i.e., daily 
RXTE observation) and very long data stretches that in- 
clude all data that fall within the decoherence timescale 
(see Section 3). The choice of different timescales is made 
to inspect the presence of rapid episodes of high ampli- 
tude pulsations (which would be missed in long-time av- 
erages) or whether pulsations are continuously present 
but with a very low amplitude (by averaging large data 
stretches) . The first folding iteration uses the orbital and 
spin solution reported in Kaaret et al. (2006). We then fit 
our pulse profiles with a sinusoid plus a constant to deter- 
mine the pulse time of arrivals (TOAs) and their pulsed 
fractions. We set the confidence level for the detection 
of pulsations at 3.6tr defined as the ratio between the 
pulse amplitude and its statistical error (see Patruno et 
al. 2010). This value is chosen to guarantee less than one 
false pulse detection when considering the entire amount 
of trials {N t „ a i s ps 3000). We then fit the TOAs detected 
until MJD 53,602 with a Keplerian circular orbit and a 
constant spin frequency v s with the software TEMP02 
(Hobbs et al. 2006; Edwards et al. 2006) and repeat the 
entire folding procedure until we reach convergence for 
our timing solution. The reason why we fit the TOAs 
until MJD 53,602 is that the spacing between detected 
pulses is sufficiently dense to avoid over-fitting of the or- 
bital parameters (see for example Hartman et al. 2009 
for a discussion). 

3. RESULTS 

The orbital solution of HETE J1900. 1-2455 is reported 
in Table 1. We find consistency between our results and 
those of Kaaret et al. (2006) within the statistical un- 
certainties, although our orbital period has an order of 
magnitude higher precision thanks to the longer baseline 
(Kaaret et al. 2006 solution refers to data up to MJD 
53,559). The precision crp OTb of our orbital period P or b, is 
such that N or b<Jp orb < P or b, where N or t, ~ 40,000 is the 
number of orbits that HETE J1900. 1-2455 completes in 
7 years. Therefore we can extend our orbital solution to 
the entire baseline spanned by the RXTE observations. 

In principle, v s is also known with sufficient precision 
from Kaaret et al. (2006) that we can confidently predict 
the pulse phases at any given epoch. However, there is 
a fundamental complication in this case represented by 
the poor knowledge of the spin frequency derivative P s . 
The presence of a spin frequency derivative has a par- 
ticularly dominant effect in HETE J1900. 1-2455 given 
its long observational baseline. We also need to consider 
the presence of timing noise and the systematic errors 
associated with the X-ray position of the source as given 
by Chandra Fox (2005). This introduces a spurious pulse 



TABLE 1 

Orbital Solution of HETE J 1900. 1-2455 



Orbital period, P orb (s) 4995.2630(5) 

Orbital period derivative, P ovb (10~ 10 s s -1 ) a <1.2 

Projected semi-major axis, a x sini (light-ms) 18.44(2) 

Time of ascending node, T asc (MJD, TDB) 53549.130943(9) 

Eccentricity, e (95% confidence level) < 4 X 10 — 3 

frequency derivative v p with magnitude of approximately 
— 10~ 14 Hzs _1 that needs to be taken into account. The 
variability of v v induces decoherence of the signal which 
is proportional to the strength of the pulse frequency 
derivative i>„: 



This value gives a maximum baseline of w 40 d and « 120 
d for v p = 10 _13 Hzs _1 and v v = 10 _14 Hzs _1 , respec- 
tively. Therefore we cannot fold very long data stretches 
without smearing the pulsations (if present). Further- 
more, even assuming that v s is zero and no timing noise 
is present in the data, we are limited by the decoherence 
time introduced by the spurious pulse frequency deriva- 
tive due to the source positional error, which is of the 
order of 100 days. 

3.1. New Pulse Episodes 

We detect 4 new pulse episodes after MJD 53,602 with 
an amplitude between 1% and 0.5% rms and one between 
MJD 53,584 and MJD 53,596 with an amplitude of 0.3% 
rms. We also find two marginal detections between MJD 
54,865 and 54,967 (3cr) and between 55,081 and 55,170 
(3.2cr). The two marginal detections have fractional am- 
plitudes of 0.1% and 0.15% rms, respectively. However, 
we do not include them in our analysis since they need 
confirmation. Robust detections are made in segments 
of different length, with higher amplitudes found in spo- 
radic and short data segments whereas the lowest am- 
plitudes are found in observations-long data segments. 
The new pulse episodes are not detected in coincidence 
or close to bursts observed by RXTE, with a minimum 
time interval of a few days between the pulse episode 
and its closest burst. The first pulse episode after MJD 
53,602 is detected at MJD 53,624, whereas the last one 
detected appears at MJD 54,499 giving a baseline of more 
than 2.5 years to perform coherent timing studies. Upper 
limits on the pulsed fractions range from 1% rms down 
to 0.05% rms (95% confidence level) with most of the 
segments having upper limits of about 0.5% rms. 

3.2. Spin Frequency Derivative 

All AMXPs discovered so far show erratic variations 
of the pulse phases on timescales of a few hundred of 
seconds up to several months that seem at first sight 
unrelated with true spin variations of the NS. These 
phase variations, commonly called X-ray timing noise, 
have been found to be correlated with variations of the 
X-ray flux in at least six AMXPs (Patruno et al. 2009b). 

To verify whether such a correlation is also present 
in HETE J1900. 1-2455 we use the correlation coherent 
analysis described in Patruno (2010b). With this method 
we fit a pulse frequency and frequency derivative to the 
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Fig. 1. — Left Panel: Phase residuals of HETE J1900. 1-2455 calculated up to MJD 53,582 when pulsations are sufficiently densely 
sampled. Each point is an observation long average. The residuals arc calculated with the correlation coherent technique (see main text). 
Central Panel: 2-16 keV X-ray lightcurve of the same data. Note that the plot shows only data in which pulsations are detected. Right 
Panel: Anti-correlation between the pulse phase residuals and the X-ray flux. 



TOAs and then choose the two parameters that mini- 
mizes the x 2 °f the linear fit between phase and flux. 
This is different than choosing the v s and i> s that min- 
imize the pulse phase residuals as it is usually done in 
standard coherent timing analysis. We find that there 
is a clear anti-correlation between the pulse phase and 
the X-ray flux. The data points exhibit a tight anti- 
correlation when considering data points up to MJD 
53,582 (see Figure 1). The x 2 of the fit although still 
unacceptable (x 2 /dof w 4, with 23 dof) is much bet- 
ter than what can be obtained with standard coherent 
analysis (x 2 /dof ~ 11 with 23 dof). When considering 
the whole data collection (including the sparse detections 
up to MJD 54,499) the anti-correlation is still present 
but with larger scattering and a worse overall x 2 (but 
still significantly better than what can be obtained with 
standard coherent analysis). 

The reason why the anti-correlation becomes worse af- 
ter MJD 53,582 might indicate that v s is not constant 
as we are assuming in our fit. We therefore split the 
data in seven overlapping intervals of different length 
and measure the v s and v s in each segment with the 
correlation coherent analysis method. We are forced to 
use overlapping intervals because the data quality is not 
sufficient to allow a measurement of v s for independent 
non-overlapping intervals. Before fitting each interval we 
change the reference epoch of our ephemeris so that each 
v s and i> s refers to a different epoch that is representative 
of the interval we are fitting. The errors are calculated 
by multiplying by \J x 2 /dof the statistical errors corre- 
sponding to the 90% confidence interval. 

We find that v s is clearly changing over time and that 
it is well described (x 2 /dof = 5.9/4) by fitting an expo- 
nential decay law with a constant level baseline: 

v s {t) = AExp(-t/r) + C (2) 

with an e-folding time r = 26 ± 4 days and a constant 
level C corresponding to a v s = 4.2 x 10 -13 Hzs -1 . The 
behavior of v s (t) is well described by the integral of the 
spin frequency derivative: 

v s (t) = J v s (t) dt = -ArExp (-t/r) + Ct + K (3) 

where K is a constant of integration that represents v a 
at the beginning of the outburst. The results for v s are 
shown in Figure 2. Given that v s changes over time we 
re-fold our data and refit the TOAs considering the dif- 
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Fig. 2. — Evolution of the spin frequency derivative with time. 
The spin frequency behavior is consistent with an exponential de- 
cay with an e-folding time r = 26 ± 4 days and a base level 
C = 4.2 X lO-^Hzs" 1 . 

ferent strength of v s in different intervals but we find no 
additional pulse episodes. 

3.3. Enigmatic Pulse Phase Drift 

A very interesting feature is evident in the timing resid- 
uals in coincidence with the flare at MJD 53,559. In 
observation-long pulse profile averages, the pulse phase 
appears offset by about 0.3 cycles with respect to the 
phases before the flare (see left panel of Figure 1). After 
the flare the pulsations disappear for several days and 
reappear aligned with the phases before the flare. When 
inspecting 300s to 500s-long pulse profile averages, the 
pulse phases during the flare show a very fast evolution. 
The phase starts 0.7 cycles (250°) offset with respect to 
the average pulse phase measured in the previous ob- 
servation (which corresponds to phase in Figure 3) and 
then drifts back by about 0.5-0.6 cycles (~ 180°) towards 
the phases of the pre-flare observations. In Figure 3 a 
negative/positive phase residual means that pulsations 
arrive earlier/later than predicted by the model. An off- 
set of -0.7 cycles is equivalent to an offset of +0.3 cycles, 
but the interpretation with the negative sign is the cor- 
rect one because the pulsations show a linear drift with 
each successive pulsation lagging the previous one. The 
timcscale for the drift is « 3000 s and pulsations have 
fractional rms amplitude which remains approximately 
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Fig. 3. — Pulse phase drift observed at MJD 53,559. Each data 
point corresponds to folded lightcurve segments of 300-500 sec- 
onds for a total observation time of about 3000 seconds. Nega- 
tive/positive values correspond to pulsations arriving earlier /later 
than predicted by the timing model. 

constant within the statistical errors with a slight excess 
in the middle of the observation. There, the fractional 
amplitude reaches 1.6% rms and stays at 1.1% rms in the 
rest of the observation. The pulsations disappear after 
this event until MJD 53,573 with 95% c.l. upper limits 
between 0.35% rms and 1.5% rms. The relatively larger 
error-bar of the pulse phase obtained by folding the en- 
tire data record at the high flux point (as shown in Figure 
1) is therefore artificial because the pulse amplitudes are 
partially smeared out by the fast pulse phase drift. 

4. DISCUSSION 

The detection of sporadic pulse episodes in the 2.5 
years following the last robust detection at MJD 53,602 
strengthens the suggestions of Galloway et al. (2008) that 
pulsations might be always present in HETE J1900.1- 
2455. The pulses we detect might represent the "tip of 
the iceberg" of very weak pulsations that are present at a 
level of < 0.1% rms, since our upper limits reach the most 
stringent value of about 0.05-0.5% rms. This suggestion 
is reinforced by our detection of pulses with amplitudes 
as low as 0.3% rms. Such low values can be reached only 
because we have an initial timing solution for the orbit 
and the NS spin which is sufficiently precise to allow a 
coherent analysis over the entire 7 years of observations. 
It is not possible, with current instrumentation, to in- 
spect any other non-pulsating LMXB with a comparable 
X-ray flux and place such extreme upper limits for the 
pulsed fractions with incoherent timing techniques. Ex- 
isting upper limits on the pulse fractional amplitudes in 
non-pulsating LMXBs arc usually of the order of 1% rms. 
The baseline available for measuring the spin frequency 
and the spin-up is ~ 2.5 years and this is an unprece- 
dented possibility for AMXPs, since accretion torques 
have been investigated so far only on relatively short 
timescales reaching at most ~ 100 days (Patruno et al. 
2010). The spin up of HETE J1900. 1-2455 follows an 
exponential decay with an e-folding timescale of 26 ± 4 
days and a constant baseline of about 4 x 10 -13 Hzs -1 . 
The spin frequency derivative expected from accretion 
torques that develop at the disk-magnctosphcre bound- 
ary of standard thin disks truncated by a constant dipo- 



FlG. 4. — Evolution of the magnetic field according to Eq. (4) for 
M_io = 10~ 9 Mq yr — 1 and an exponentially decaying i> a as shown 
in Figure 2. The vertical dashed blue line represents the time when 
pulsations first disappear and corresponds to a field of 7 X 10 7 G. 
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where M_io is the mass accretion rate in units of 
10 Mo/yr, Bs the dipolar magnetic field at the poles 
in units of 10 8 G, Mi. 4 the mass of the NS in units of 
1.4M©, Rio its radius in units of 10 km. The param- 
eter £ is introduced to account for the uncertainties in 
the torque at the edge of the accretion disc and is in the 
range £ « 0.3 - 1 (Psaltis & Chakrabarty 1999). 

As noticed by Galloway et al. (2008), the lightcurve 
of HETE ,11900.1-2455 shows an erratic behavior with 
large variations in X-ray flux on timescales of tens of 
days. The flux averaged over several tens of days is, how- 
ever, rather constant (we find variations of at most 20% 
in X-ray flux between the averages of our seven inter- 
vals) and is certainly not exhibiting an exponential drop. 
Therefore the exponential decay of i> a cannot be related 
solely with variations in M (if we assume that Lx oc M; 
see, however, van der Klis 2001 for a criticism), neither 
with variations of M and R which stay basically constant 
throughout the outburst. If the magnetic field is respon- 
sible for the variation of v Sl then the B field has to decay 
approximately exponentially with a very short e-folding 
time of ~ 10 days. 

Magnetic field burial models predict an exponential de- 
cay of the external magnetic field which is screened by 
freshly accreted plasma (Cumming et al. 2001; Cumming 
2008). As matter is accreted to the polar caps it will 
eventually spread laterally and bury the field underneath. 
If the new material accumulates on the NS surface on a 
timescale which is much shorter than the Ohmic diffusion 
timescale, the magnetic field is (partially) screened. In 
these models the field suppression operates on timescales 
that vary significantly and depends on several assump- 
tions. For example, Choudhuri & Konar (2002) have 
shown that screening timescales of ~ 1 yr can be achieved 
provided that the no magnetic buoyancy is present (see 
also Payne & Melatos 2004, 2007; Wette et al. 2010). 
The fact that our fit requires a constant level can also 
be naturally explained in this scenario by considering 
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that the magnetic field is not efficiently screened once its 
value is so low that channeled accretion becomes diffi- 
cult (see Konar & Choudhuri 2004 and compare our Fig- 
ure 2 with their Figure 10; Zhang & Kojima 2006). Once 
the outburst is over, the magnetic field can re-emerge 
on the Ohmic diffusion timcscalc (Cumming et al. 2001) 
and return to its initial value of ~ 10 8 G. For illustra- 
tive purposes we plot in Figure 4 the B field at the NS 
poles of HETE J1900. 1-2455 where we use Eq.(4) and 
assume a constant M_io = 10, £ = 1, M\a = 1 and 
i?io = 1. The AMXP has a magnetic field at the begin- 
ning of the outburst of B L ~ 5 x 10 s G and a final field 
of Bf ~ 7 x 10 5 G, which is significantly less than the 
minimum field necessary to truncate the accretion disk 
(Psaltis & Chakrabarty 1999). This could in principle 
be related to the extraordinary low rms amplitudes of 
pulsations in HETE J1900. 1-2455. 

A similar behavior might not be observed in other 
persistent AMXPs possibly because of the substantially 
smaller mass accretion rate which is one to two orders 
of magnitude lower than in HETE J1900.1-2455. If the 
timcscale for the magnetic screening scales inversely with 
the mass accretion rate (Konar & Choudhuri 2004) then 
the magnetic field of persistent AMXPs might require 
one to several years to substantially decrease. Since 
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the outburst duration of persistent AMXPs is at most 
100 days, the screening mechanism cannot reduce the 
strength of the magnetic field below the level necessary 
to channel plasma along the field lines. 

If magnetic screening is the correct interpretation of 
our findings, then the anti-correlation between phase and 
X-ray flux might possibly also be explained with varia- 
tions of the magnetosphere (and thus on the position of 
the polar caps) that responds to variations in the amount 
of accreted material (see also theoretical investigations 
of this problem in Long et al. 2012; Kajava et al. 2011; 
Poutanen et al. 2009; Lamb et al. 2009a). However, the 
anti-correlation does not work on very short timcscalcs 
of the order of a few hundred seconds. In particular, 
the sudden w 180° (0.5 cycles) pulse phase drift ob- 
served at MJD 53,559 remains an enigmatic event and 
it might contain the key to understand why pulsations 
disappeared for the first time right after this event. 
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